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PACS numbers: 



We propose a new model of slow-roll inflation in string cosmology, based on warped throat super- 
gravity solutions displaying 'walking' dynamics, i. e. the coupling constant of the dual gauge theory 
slowly varies over a range of energy scales. The features of the throat geometry are sourced by a rich 
field content, given by the dilaton and RR and NS fluxes. By considering the motion of a D3-brane 
probe in this geometry, we are able to analytically calculate the brane potential in a physically inter- 
esting regime. This potential has an inflection point: in its proximity we realize a model of inflation 
lasting sixty e-foldings, and whose robust predictions are in agreement with current observations. 
^ I We are also able to interpret some of the most interesting aspects of this scenario in terms of the 

U , properties of the QFT dual theory. 

(N 

>: 

I. INTRODUCTION 

■ The inflationary paradigm provides convincing solutions to the basic problems of standard big bang cosmology. 
„ Moreover, it offers a testable mechanism for the generation of primordial cosmological perturbations that seed the 
' formation of large scale structures and are imprinted on the cosmological microwave background radiation [1]. So 
far, the simplest models of inflation are in agreement with observations; they are based on a single scalar fleld, the 
inflaton, which slowly rolls down a flat potential for a period sufliciently long to achieve sixty e-folds of inflation. 

• ■ 

i While at flrst sight it seems easy to build inflationary scenarios with the desired properties (and indeed there are 
^ ' plenty of proposals), a theoretical challenge is to embed consistent models of inflation in a fundamental theory, capable 
5^ , to combine cosmological inflation with particle physics models at short distances. An option is to embed inflation in 

string theory, using as inflaton candidate one of the many light moduli available in string set-ups. See [2] for recent 

nice reviews that also refer to the relevant literature. 

Much work has been done in particular for embedding models of inflation into warped flux compactifications of type 
IIB string theory. This approach is referred to as warped D-brane inflation [3], in which inflation is realized by the 
motion of a probe D3-brane along the radial direction of a strongly-warped throat region inside a flux compactiflcation. 
In this type of geometries, fluxes [4] and non-perturbative effects, such as gaugino condensation on wrapped D7-branes 
[5], are turned on to stabilize light moduli that do not drive inflation. This approach has been reflned significantly 
over past years, using the deep understanding of warped conifold geometries recently acquired, motivated by the 
AdS/CFT correspondence. 

In this work, we continue to build along this line of research, embedding a slow-roll model of inflation in a warped 
throat geometry with various background flelds and fluxes turned on. We will focus on a warped throat dual to a 
gauge theory that exhibits 'walking behavior' [6, 7], i. e. which has a coupling constant characterized by an RG flow 
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that runs slowly within a range of energy scales. From a geometrical point of view, this energy range corresponds to 
a region of transition between a geometry sourced by a stack of D5-brancs in the deep infrared (IR), and a space-time 
that approximates the Klebanov-Strassler (KS) throat [8] in the ultraviolet (UV). This region is controlled by fluxes 
that acquire non-trivial profiles which are tunable to a certain extent. The rich content of background fields turned on 
in the geometry , which is necessary to generate multiple scale dynamics in the QFT dual, offers an attractive arena 
for embedding models of warped D-brane inflation. This is the scope of this work, in which we show that a successful 
inflationary model can be embedded in this geometry using a probe D3-brane. Fluxes and warping generate a force 
which acts on the probe D3-brane, thereby inducing inflation with no need to include additional sources such as anti- 
D3-brancs: in this sense, our scenario can be seen as a generalization of the inflationary model built on the baryonic 
branch of the Klebanov-Strassler throat [9]. In our set-up, slow- roll inflation naturally occurs precisely in the region in 
which the walking dynamics is manifest in the dual fleld theory; hence the name slow-walking inflation. Remarkably, 
in the regime we are focussing on, we find a fully analytical form for the probe brane potential, which enables us to 
clearly appreciate the roles of the geometrical parameters that govern this potential. From the perspective of the dual 
field theory, it is possible to identify which operators are responsible for generating the main features of the potential, 
leading to a QFT interpretation of some aspects of our scenario. In particular, the VEVs of two operators in the 
QFT dual start to dominate in the walking regime, and combine in such a way to generate the flat potential for the 
D3-brane. The large number of parameters characterizing the geometry provide enough freedom to tune the D-brane 
potential, allowing to comfortably accommodate the requirements of producing sixty e-folds of slow-roll inflation, and 
to match the observed amplitude and scale-dependence for the power spectrum of scalar fluctuations. Interestingly, 
inflation occurs naturally around an inflection point of the brane potential, which can be related to the non-trivial 
proflle of the dilaton along the radial direction. Within our scenario, we obtain robust predictions for the properties of 
the inflationary process: the slow-roll e-parameter is much smaller than the ry-parameter; the Tj-parameter is negative, 
and its absolute value is proportional to the inverse of the number of e-folds. The resulting value for the scalar spectral 
index is in good agreement with observations by WMAP [10]. 

Building a model of inflation in the warped throat geometry, as we achieve here, is only the flrst step in the theoretical 
investigation of our scenario. Next, it needs to be checked whether inflation survives after including corrections to the 
D-brane potential due to couplings between the D3-brane, the moduli stabilizing sectors and gravitational degrees 
of freedom as well as Kaluza-Klein modes. This is the so-called r?-problem: typically, corrections to the D3-brane 
potential due to the aforementioned effects tend to make the 77-parameter large, and spoil inflation. Originally, this 
problem was pointed out in [11] when attempting to embed inflation in supergravity. It has been widely analyzed 
in recent years in the context of warped D-brane inflation [3, 12-23], when embedding a D3-brane in the AdS-like 
region of a Klebanov-Strassler throat. Our set-up is different from the KS warped throat D-brane scenario. Thus, a 
careful analysis of the Tj-problem in the present context is not straightforward so we will not pursue it here. Let us 
however point out that the large number of parameters involved and the various operators that can be switched on 
from the QFT perspective, leave the hope that dangerous contributions from compactiflcation effects can be tuned 
away without qualitatively changing the features of our set-up. 

The plan of the paper is the following. In Section II, we carefully describe the geometry that constitutes the arena for 
our inflationary model, with emphasis on the parameters describing the warped throat (see in particular Section II D). 
In Section III, we briefly review D3-brane inflation in a generic geometry of the above form. In Section IV, these results 
are then applied for the background of Section II. We find (for a particular range of parameters) analytic expressions 
for the quantities for inflation (e, rj, and the amplitude of the power spectrum), showing that they exhibit a scaling 
behavior with the parameters of the model. An explicit inflationary trajectory is also presented. Finally, Section V 
summarizes our results and concludes pointing to future possible studies. Various appendices of technical nature are 
included. We note that other models of inflation based on embedding probe branes within gauge/gravity duality have 
recently been discussed in [24-30] . 

II. THE WALKING GEOMETRY 

In this section, we discuss the warped throat geometry of interest, emphasizing its geometrical and physical aspects 
that we will use later for building our model of inflation. The same class of geometries were considered in [6, 7] , where 
they were instead used for building particle physics models of walking technicolor. 

A. Wrapped-D5 system 

In this section, we will consider type IIB string theory in the Supergravity limit. The topology of the background 
space time will contain a four dimensional Minkowski space and a six-dimensional Calabi-Yau (CY) 'internal' space 
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related to the conifold. Due to this, the background space-time will preserve minimal {M = 1 supersymmetry in four 
dimensions. 

Let us start then, by considering typo IIB supcrgravity and the geometry produced by stacking on top of each 
other Nc D5-branes wrapping an inside a 6d Calabi-Yau (CY) cone [31, 32] — e.g. the conifold with topology 
M. X X S^. We truncate type IIB supcrgravity to include only the metric, the dilaton $ and the RR three-form 
F3, which we express in terms of the {SU{2) left-invariant) one-forms uii, {i = 1,2,3) sec Appendix A for more 
details. We then propose an ansatz that assumes the functions appearing in the background to depend only on the 
radial coordinate p. Setting a'gs = 1, we write the background metric in Einstein frame as 

H — [(wi -I- a{p)def + {Cb2 - a{p) sm.ed4>f] H ^(^3 + cos6»#)^ 

As mentioned above, there are also a dilaton $ and a RR-three form F3 that complete the background. We write 
details of the full configuration in Appendix A, which as we commented, will preserve minimal SUSY. 

Any particular SUSY solution is determined by solving the system of BPS equations. These equations are in general 
non-linear and coupled. They can be rearranged in a convenient form by rewriting the functions appearing in Eq. (1) 
in terms of a new set of functions (a 'basis' that we call P, Q, .... defined in Eq. (A. 3) in Appendix A), to make the 
system still non-linear but at least decoupled. 

After some lengthy algebra (described in detail in [32] and summarized in Appendix A), we find that the relevant 
functions for the purposes of this paper are given by 

g(p) = iV,(2pcoth(2p)-l), 2e"^ = P\ e^^ = ^-^^^^^, (2) 

with $0 an integration constant and the function P{p) satisfying the second order differential equation 

We will refer to Eq. (3) as the m,aster equation. This is the only equation that needs solving in order to generate 
the large classes of supersymmetric background solutions we are interested in. The equation for P depends on the 
function Q defined in Eq. (2), proportional to the number Nc of wrapped branes. 



B. Rotation: U-duality as a solution-generating technique 

We now take the configuration of Eq. (1) and apply to it a solution generating technique. This will produce a more 
interesting background in view of our applications (more on this in subsection II D). In [33], the authors proposed 
a U-duality that takes a background of the form written in Eq. (1) with certain particular solutions of Eq. (3) and 
maps them into another background where new fluxes are turned on. See also [34] for a different perspective on this 
technique. 

The effect of this solution-generating technique (which we will refer to as 'rotation') can be summarized by defining 
a basis of 1-forms, which in Einstein frame reads 

e^' = /i ^e^dxi, = klh^e^^^dp , e'^ = k^h^ — — {uj^ + cos 9d(p) , 

= klhiei+''de, e"^ = klhiei+'' smOdip , (4) 
e =klh* — ^ — (Cji + adO), e =klh^ — ^ — {Cj2 — asmOdip) , 

where Xi arc the four Minkowski directions (and the uJi are defined in Eq. (A.l)). We define the warp factor h 
appearing in the previous expression as 

h^k^^l-kle^^) , (5) 

which depends on two constants, ki and ^2. These will be determined later on by appropriate physical requirements. 
The new configuration contains a metric Qmn, a dilaton a RR three-form and a five- form F5 together with a 
NS three- form and its associated potential B2. 
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The metric and the RR five form — the two most relevant quantites for our purposes — are written in terms of 
the vielbein of Eq. (4) as (with e*-' - '' = e' A e-* A • • • A e'') 

10 

ds| = ^{e'f = h'^e^dxl-i + kih^e^+'^^dp^ + ... , 

(6) 

dp \ h J 

The main message up to this point is the folfowig. We start with a configuration of D5-branes preserving minimal 
SUSY. This set-up is characterized by the solutions to what we called the master equation, Eq. (3). Some solutions 
to this equation can be used to generate other backgrounds like the ones in Eq. (6) — with full details given in 
Appendix A. Regardless of the background being that in Eq. (1) or that in Eq. (6), there is one master-equation 
controlling them. Once we have solved for P{p), we can study the dynamics of probe D3-branes embedded in the 
backgrounds of Eqs. (1) or (6) — passing from one to the other is just a matter of appropriately choosing the constants 
fci, k2. Therefore, we will find useful to present a classification of the solutions to Eq. (3) in the next section. 



C. Solving the master equation: walking solutions 

We now characterize the properties of the solutions of Eq. (3). Most of the results of this Section were already 
discussed in [6, 32], hence we will be sketchy. 

Since the master equation is a second order differential equation, the general solution to Eq. (3) will have two 
integration constants. In a regime in which P ^ Q (we will discuss its physical meaning below), the master equation 
is approximately solved by 

~ 1 /3 

P = c [cos'^ a + sin^ a (sinh(4p) - 4p)] , (7) 

which depends on the two integration constants a and c. By inspecting the previous expression, one notices that it 
is possible to determine a scale p* such that the solution for P is approximately constant for p < p^, while for p > p^, 
behaves as P oc e~ . This intermediate scale p* will play an interesting role below. In the limit of small a, one finds 
that approximately [6] 

4p, ~log(2cot^a) . (8) 

This allows us to rewrite 

r . ii/3 
P ~ csina e^P' + 2(sinh(4p) - 4p) . (9) 

Starting from this observation, a more precise characterization of the solution for P proceeds as follows. 
In the far UV, for p — > oo, one can check that the following expression solves Eq. (3) 

P = c^e-ol^ + 4^ (p^ - p + + {-\c^p - i^^) e'^ol^ + 0(e-^) , (10) 

where c± are the two constants characterizing these solutions. While it will be clearer below, notice that by comparing 
this with the large-p limti of Eq. (9) we can identify c+ c sin a. 
In the deep IR, for p ^ 0, one finds 

P = CO + ksCop' + %C0P' - klcop' + 16(2c^fe3 -5fc3A^e^) ^7 + 0(^8) (11) 

iUoco 

where now cq and arc the free parameters. 

One can hence write all of these solutions by specifying Nc and for instance any of the pairs {c,a), (c-|_,c_) or 
(co, ^3)5 imposing that the function P smoothly connects between IR and UV. Which parameterization to use is 
mostly a matter of convenience. We found that for analytic computations a particularly convenient choice is (c+, p*), 
while for numerically solving the master equation we instead specify (co,p*). We show some numerical solutions to 
the master equation for the function P(p) as well as the function Q{p) in Fig. 1 and the corresponding background 
functions e^^^p\ e^'"'^"^ and h{p) in Fig. 2. Note that the warp factor h{p) has significant warping only until p ^ p*. 
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FIG. 1: Plots of P and Q as a function of p for different values of co, p* and Nc- The color coding in the left figure is as follows: 
curves with the same color have a common value of co between 100 and 1000 but differ in their value of p* ~ 1, 2 or 3, while 
A^c = 1 is always kept fixed. For the right figure, each color corresponds to a different value of A''^ from 1 to 6. 



To continue further, we need some relation between the UV and IR expansion coefficients. In [35], some approximate 
relations between the coefficients c±,cq, of the UV and IR expansion and c, a of the approximate solution in Eq. (7) 
were derived. We can use them to express Co, and c_ in terms of c+ and p* (at leading order in an expansion in 
large c+ and p*): 




FIG. 2: Plots of e"*, e^^ and h as a function of p for different values of Co and p* with Nc fixed. The color coding is as follows. 
In the left and middle figure, curves with the same color have a common value of co between 100 and 1000 but differ in their 
value of p, ~ 1, 2 or 3, while TVc = 1 is always kept fixed. In the right figure, curves with the same color have the a common 
value of p, ~ 1 (blue), 2 (green) or 3 (red) but have a different value of puv = 10 or 12 at which we impose the normalization 
condition /i(puv) = 1, see Section II D. 



More systematically, we can describe a solution to the master equation (3) using a perturbative approach. The 
solution is then formally given in terms of an infinite series of non-explicit integrals. The usefulness of this formal 
method will become clear in what follows. Since this strategy for solving the equation was developed in [6] — see 
also the Appendix B of [32] — we will just summarize here the main points. Note that this approach is essentially a 
systematic expansion around the (approximate) solution in Eq. (9). 

We can integrate twice the master equation (3), taking already into account the UV and IR asymptotics described 
above, to get 

/•p rp roc, P'OO' f 

P^-3Q^P + 6 dpPQQ'-U dpsmh^{2p) / dp f-^ = 16ci^/ dp sinh^ {2 p) + P{Of , (13) 

Jo Jo Jp smh (2p) Jo 

The double integration provides two integration constants, P(0) and c+. The latter is related to the UV behaviour of 
P(p) we discussed above, P{p) — c^e^P^^ + . . . . Whilst P{Q) is related to the IR behaviour of P{p)- For the walking 
solutions considered here, we have P{Q) = cq = c+ e^'''!'^ . 
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We can formally solve Eq. (13) in a series expansion in c , ^. We propose the ansatz 



Pip) 



n=-l 



'P-n = C+Pi + Po 



-P-1 



-P-7 



Inserting this ansatz into Eq. (13) and matching order by order in c+ one finds that Pq — P-2 — 
P, = (2 (sinh(4p)-4p) + , 



52 dpPiQQ'-A dpsinh^2p dp - 

1 V Jq Jo J 5 si 



PiQQ' 



sinh 2p 



(14) 

= P-2k = and, 
(15) 
(16) 



where we have used P(0)^ 



^3 g4p, infiicated in Eq. (12). Note again that in Eq. (9) we can identify 



c+ ~ c sin a. Note also that up to the factor c sin a ~ c_|_ Pi is nothing but P from Eq. (9). A recurrence relation for 
P_(2fe+i) can be found by setting TV/ = in Eq. (B.7) of [32]. This series converges rapidly to the numerical solutions 
of Eq. (3) — see Fig. 3. However, the information contained in Eqs. (15) and (16) is already sufficient for our needs. 





FIG. 3: Plots of P and P' as a function of p for Nc = 1, co = 100 and p, ~ 2. We show the result of solving the master equation 
numerically (red dots) as well as three approximations using either only the leading term of the c^^ expansion P ~ c+Pi (black), 
the leading term of the UV expansion P r . p4p/3 
(purple dashed). 



c+e ''' (blue dotted) or the leading piece of the IR expansion P ~ co + co fea p 



Before we close this section, let us briefly elaborate on the physical meaning of the approximate solution in the 
regime P ^ Q, which is given by P in Eq. (7). This form of the solution corresponds to focusing on the limit of 
large c+, neglecting the terms weighted by inverse powers of this quantity in the expansion Eq. (14), i. e. dropping all 
but the first term. It was shown in [34] that in such a limit the geometry of the internal space becomes close to the 
deformed conifold (or variations of it once we turn on the constant e''* '-^ cos^ a) . This scaling comes together with 
the dilaton approaching a constant value throughout the entire throat. In other words, for large but finite values of 
c+ we have a background which is 'close' to KS but not exactly equal to it. As we will see, this leads to a situation 
in which the force on a probe D3-brane vanishes as c_(- oo since the latter does not break SUSY in pure KS (see 
also Eq. (15.3) of [9]). 

Let us also briefiy comment on the singularity structure of the walking solutions. The backgrounds where e^^* > 0, 
do indeed contain a singularity at p = 0. This singularity is 'mild' in the sense that the Ricci scalar R, the scalar 
RabR"''' are finite, but the Kretschmann scalar RabcdR"'^'^''' diverges. Nevertheless, this singularity does not affect our 
scenario or the results derived from it. As we will discuss below, the inflationary dynamics take place in the interval 
[pi, p,], far away from the singularity. Hence, our treatment is free from pathologies and we are calculating in a set-up 
which is reasonable within the supergravity approximation. 

The reader may have observed that in this presentation of the solutions to the master equation (3) and its rotation 
to the configuration in Eq. (6), various integration constants have appeared: [fti, A;2, $oj c+, p*] are the ones relevant 
for the topics of this paper. Other constants — like the choice of origin for the radial coordinate or a constant 
appearing in the function Q{p) of Eq. (2) — either have uninteresting physical meanings, or have been fixed to values 
that avoid nasty singularities in the geometry. See [35] for a complete treatment of all integration constants. It is thus 
useful to characterize the solutions in terms of the relevant constants, imposing some physical criteria to fix some of 
them. We now turn to discuss this topic. 
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D. Physical characterization of the geometry 

The warped throat geometries described above have various features that make them interesting for appUcations to 
particle physics model building and, in the present context, for cosmology. Under certain situations, the integration 
constants characterizing the geometry have a counterpart in a dual QFT, where they can be interpreted as VEVs or 
couplings of operators perturbing a conformal fixed point. Here, we briefly review the most relevant properties for 
our purposes, referring the reader to [6, 35] for more details. 

If the integration constant p* is sufficiently large and positive, there is an intermediate regime in the radial direction, 
from /J/ ^ 1 to p ~ p», in which the effective 4d gauge coupling of the dual field theory is finite and approximately 
constant [6]. Hence, the geometry associated with solutions with large p* are suitable for describing field theories 
exhibiting walking dynamics. In the QFT dual, the VEVs of two operators start to dominate in the regime pi < p < p*, 
which arc absent in the KS configuration: a dimension-2 VEV that brings the background on the baryonic branch [36], 
and a dimension-6 VEV associated with the walking dynamics. The combined action of the two operators, enriches 
the dynamics generating the intermediate walking region and provides yet another physical scale, at the value p/ ~ 1 
in the radial coordinate p, corresponding approximately to the scale where the walking regime ends and the system 
enters into a confining regime. At the geometrical level, the position pj roughly corresponds to the value of the radial 
coordinate below which the functions a(p), b{p) (controlling F3, H3) defined in Appendix A become non-trivial — this 
fact is associated with the spontaneous breaking of a discrete Z27v^ global symmetry and with a confining behaviour 
— see the second paper of [7]. As we will see in what follows, the region p/ ^ p ^ p* is important for building our 
inflationary set-up. 

In order to consider inflationary models, we need to glue our throat geometry to a compact space, say a (warped) 
CY-three-fold, in the UV. To this extent, let us discuss the UV behavior of h{p), which is controlled by the constant fc2. 
For generic values of the parameter ^2, the warp factor h{p) asymptotes to a constant as the radial coordinate p grows. 
In the field theory dual, this is associated with the presence of a dimension-8 operator that dominates the dynamics 
in the UV and needs a careful specification of a UV completion — see the discussion in [35]. In this UV region, the 
geometry generically differs considerably from a space-time dual to a conformal theory. Thus, it makes little sense to 
interpret the duals of the integration constants of the geometry as operators representing 'small perturbations' around 
a conformal fixed point. Fortunately, the UV behavior can be improved: tuning the parameter k2 appearing in the 
warp factor h, cf. Eq. (5), to be ^2 = e~*^°°^ adiabatically switches off the aforementioned dimension-8 operator, 
and the geometry becomes (logarithmically) close to AdS for large values of the radial coordinate [35], thanks to RR 
five- form and NS three- form fluxes that get turned on at p ~ p* — see the plots in [35] . The configuration at p — > 00 
becomes almost identical to the Klebanov-Strassler background [8] . So by making this choice for /c2 , we can reliably 
identify the operators in the dual fleld theory representing 'small perturbations' from a quasi-conformal regime. This 
is good news also for the aim to build an inflationary model. After tuning k2 as explained above, in the UV, the throat 
is very similar to AdS space, which implies that upon gluing the throat into a compact space the (four-dimensional) 
graviton zero mode will be mostly localized in the compact part of the bulk and not in the warped throat. This 
information will be useful in what follows, to fix the value of the effective four-dimensional Planck mass in string 
units. Henceforth, we always set ^2 = e~*(°°^. 

Let us proceed with specifying the UV properties of our configuration. In the limit c+ 00, the dilaton becomes 
close to a constant, see for example [34, 35]. Without lack of generality, we can choose boundary conditions such 
that in the far UV $(00) = 0, such that ^2 = 1- Consider the expression for the dilaton in Eq. (2); the asymptotic 
expansion for e^* as p — >■ 00 reads 

.*-4*o ^J_f,_ ^N^e-^P/^ (8p-l) 
" 8c^ V 44 

Given the requirement above, $(00) = 0, the parameter $0 is then determined in terms of c+ as 

*o = JlnQ4) ^ e-o = ^4. (18) 

Another physically convenient requirement allows to fix the parameter ki. The warp factor h is (we set ^2 = 1), 

h = k^^{l - e^*) . (19) 

We choose the value ki such that at the UV boundary on which the throat is joined to the compact space the warp 
factor is of order one, /i(puv) ^ 1- This implies that the size of the compact bulk is determined by the string scale, 
and its volume is not enhanced by large values of the overall warp factor. This was also done in [3] and is important 
for fixing the scales that enter in the phenomenological discussion of the inflationary set-up we are going to develop. 
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Using Eqs. (17) and (18), the condition h{p\jy) = 1 implies 

/3 7V,e-4puv/3^8^yy_l 



(20) 



After making this choice, we have fixed three of the parameters controhing the warped throat ki,k2,^Q. We still can 
vary c+ (controlling the profile of the dilaton field), p*, the position puv at which the throat ends and the number 
Nc of wrapped D5-branes at the tip of the throat. 

All these parameters and scales discussed above will characterize our model of inflation. To close this section, let 
us briefly summarize the most important parameters specifying the geometry and their physical interpretation, see 
also Fig. 4. 

• c+: controls VEV of a dimension-2 operator (^^2) ^ — corresponding to motion along baryonic branch — it 
also induces dilaton proflle. 

• controls VEV of a dimension-6 operator (Og) ^ e^''* inducing walking dynamics — also modifles the dilaton 
profile, especially in IR region. 

• Puv: controls size of warp factor at the tip, i. e. ratio between confining scale and UV cutoff — distance over 
which strong warping occurs: ~ puv — P* ■ 






FIG. 4: Plots of P, e"** and h as a function of p for Nc = 1, Co = 100, ~ 3 and puv = 7. The different scales , pi and puv 
playing a role in this paper are indicated by the gray dashed lines. 



III. BRIEF REVIEW OF WARPED D3 INFLATION 



We consider a probe D3-brane placed in a warped throat geometry. The inflaton field will geometrically correspond 
to the radial position of the D3-brane on the throat (the first papers identifying the inflaton with brane motion have 
been [37-39]). To derive the canonically normalized inflationary action, we start from a generic string frame metric 
— one can show that the relevant formulae below have the same expression if we work in Einstein frame — with a 
warped conical region and RR four-form ansatz of the form 

dsl^ = Hi 3 + H2 dp'^ + . . . , C4 = Ci{p) dt A dxi A dx2 A dxs , (21) 

where the dots denote the angular part of the metric. The induced metric on a D3-brane that extends in the four 
non-compact directions and moves along the radial direction, 

1:4 = [i,xi,X2,a;3] , p = p{t) , (22) 

is given by 

dsLd = Hi (dxl + dxj +dxl) + {H2 P^-Hi) dt^ . (23) 
The BIWZ action for this probe brane is, again in string frame. 



Sbiwz — —T-i 



3 / d'^x(^e *v^--det[gj„d] - C^^j 



-Ta / d''x[e-''Ht\ll-^^p^-C, 



(24) 
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We expand for small velocities p and find 



Sbiwz = n I d'x ( ^bE2flp^ _ e-^Hl ( 1 - ) , (25) 



2 " ^ V Hi 

from which we read off the effective D3-brane potential 

V = T,e-^Hl • (26) 

The canonical radial variable, in string frame, is determined from 

dr = ^ie-^l'^ ^fU^^dp. (27) 

The determination of the potential for the D3-branc may then proceed following a two-step procedure. In a first step, 
we can find supersymmetric solutions for the background field equations, as we did in the previous section, and embed 
a probe D3-brane in the corresponding configuration. If the brane breaks SUSY, it feels a force that makes it move. 
This motion induces a modification of the warped geometry, which we assume to be limited to change the 4d metric 
dxf 3 describing the observed four dimensions, lying inside the warp factor in Eq. (1) — or in its 'rotated' version 
in Eq. (6). Namely, we assume that this 4d metric changes from Minkowski to an homogeneous and isotropic FRW 
space time. For this approach to be valid, we make the hypothesis that all the relevant moduli are stabilized with 
sufficiently high masses such that they are not destabilized by the brane motion. The second step [13-15] consists 
of starting with a fiducial warped throat configuration, for example the KS throat, on which a probe D3 feels no (or 
only a small) force, and calculating perturbatively the D3-brane potential around this configuration generated by the 
backreaction on the geometry from corrections sourced e. g. by SUSY breaking contributions in the bulk necessary to 
stabilize the geometric moduli. 

In any case, the D3-brane action then has to be supplemented by coupling it to gravity. This is done by gluing 
the throat to a compact space, which induces an Einstein-Hilbert term completing the D3-brane action. Usually, this 
gluing also provides new contributions to the D3 brane potential, whose combination is neglected in the first step 
described above, while it can be computed perturbatively in the second step at least in principle. 

The main goal of the first step, which we will pursue in this paper, is to clearly identify the forces that induce 
the brane motion, offering neat physical and geometrical insights for understanding inflation in geometrical terms, 
and possibly leading to a fully calculable inflationary potential. The results can also be interpreted in a dual QFT, 
offering interesting new perspectives on possible mechanisms that drive inflation from the; point of view of strongly- 
coupled field theories. The problem with this first step, is of course that one does not take into account all possible 
contributions to the inflationary potential since contributions depending on couplings with other light moduli or fields 
in the compact bulk are assumed to be negligible. This is in general not consistent and is generally associated with 
the so called t] problem (more on this in Section V). 

The purpose of the second step is then to offer a systematic way to calculate perturbatively the brane potential, 
offering fully reliable results, although the calculations are not always straightforward, and a further numerical analysis 
of the inflationary trajectory is often required (see e.g. [17-19] for the case of the KS throat). Ideally, the physical 
and geometrical insights about the forces acting on the brane within a given set-up, acquired by the first step, should 
be supplemented by a careful determination of the corrections to the inflationary potential, calculated following the 
second step. In this way, one can obtain a fully-controlled satisfactory scenario for warped brane inflation. 

After clarifying these methodological issues, let us specify two quantities — the slow-roll parameters e and t] — 
which are important for characterizing our inflationary model. In the limit in which the brane moves slowly, the 
slow-roll parameters are given by 



d'^V (dpV (dV\ fd?p 



Ml 
V 



dp^ \dr J \dp J \ dr^ 



(28) 



Note that these are simply the usual slow-roll parameters in terms of the canonically normalized radial coordinate r 
obtained from Eq. (27), namely 

Given a particular warped geometry and its fleld content, the resulting inflationary potential built along these lines 

can exhibit regions that are flat enough for supporting a sufficiently long period of slow-roll inflation. In this set-up, we 
can interpret the properties of the inflationary scenario using the geometrical features of the throat, or the properties 
of its dual QFT. 
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IV. INFLATIONARY CASE STUDY — WALKING SOLUTION 



Let us now discuss how inflation occurs in the geometry of the walking solution discussed below Eq. (7). In the 
string frame, the functions i?i,2, C4 in Eq. (21) for this solution are given by 



Tj 1.-1/2 * TT I t.1/2 2fe+$ /I f-l 2* ; 



(30) 



The inflationary potential V and canonically normalized radial variable r (which acts as the inflaton field) are 

V = and dr = v/JVfcTe'^+^/^dp. (31) 

e-'P + fc2 

In the following, we take A;2 = 1 and ki given by formula (20) as discussed in Section II D. 






FIG. 5: Plots of V , V' and V" as functions of p rescaled by appropriate factors. The color coding is as follows: each color 
corresponds to a different value of p* ~ 1 (red), 2 (blue) or 3 (green); in the left figure curves with the same color differ by 
their value of co = 100 or 200. We always keep Nc = 1 fixed. 



In the limit ^^/c+ — ^ 0, our geometry becomes (almost) identical to the Klebanov-Strassler throat: it would coincide 
with KS in this limit, if additionally we consider the limit — > —00; see Eq. (15). For ^^/c^ — >■ 0, the dilaton $ 
becomes a constant, and thus the brane does not feel any force. Outside this limit, a non-trivial profile for $ forces the 
brane to move, leading to cosmological evolution compatible with inflation. The main advantage of this framework is 
that a controllable brane motion is associated with the properties of our background, which is free of singularities in 
the region of interest, and we do not need to add anti-D3 objects to make the brane move. In this sense, our scenario 
can be considered as a generalization of the inflationary D-brane model in the baryonic branch of the KS throat [9] . 
Let us briefly explain the differences between the two setups. First, in [9] it was proposed to have inflation in the 
AdS-like UV region, while in our case inflation takes place near the tip in the IR. Second, as we will see below, in 
our case 77 becomes small in the vicinity of an inflection point and additionally the walking dynamics lead to a much 
stronger suppression of e by a factor of e"^''*/'^. 

Before studying in detail the quantitative features of an inflationary model that can be build using this geometry, 
let us start with some general qualitative considerations. In the deep IR, p — > 0, the geometry is essentially controlled 
by the wrapped D5-branes. The corresponding warping creates the attractive force felt by the moving D3-brane. 
Moreover, the fields sourced by the D5s, which are concentrated in the IR region, make this attractive force stronger 
and stronger when increasing the radial coordinates. On the contrary, at larger values of p, we are leaving the 
region in which the wrapped D5-branes dominate, and the warping changes. Additionally, for p > p*, new fiuxes 
(the RR five-form and NS three-form) appear and change the throat geometry, making it similar to the KS throat. 
Consequently, the attractive force acting on the brane becomes weaker, until it asymptotically vanishes in the far UV 
when the geometry quickly approaches KS. In this region, however, the throat will be attached to a compact space as 
we explained above, and drastic modifications of the background have to be expected in the gluing region. See Fig. 5 
for plots oi V, V and V" for different values of cq and p,. 

On general grounds, we expect a region of transition, in which the attractive force acting on the brane reaches a 
maximum in size at a certain value of the radial coordinate after which its size starts to decrease (see the middle 
figure of Fig. 5). In the following, our aim is to identify this region, and we will focus on it for building a model 
of slow-roll inflation with the necessary properties to match observations. As we will discuss in detail, this region is 
characterized by a potential that is flat enough to provide 60 e-folds of inflation. Remarkably, we will show that our 
set-up allows to determine and analyze the inflationary potential analytically, at least in a physically well-motivated 
limit. In the QFT dual, this region is precisely the walking region in which the coupling constant varies only mildly. 
(Recently, similar considerations have been made in a conceptually different set-up in [25].) 
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A. The D3 brane potential in the Iarge-c+ and large-p* limits 

As explained above, in the large-c+ limit, our geometry approximates the KS throat, and we expect that in this limit 
the force acting on the D3 becomes small and capable of supporting a sufficiently long period of slow-roll inflation. 
The limit of large c_|_ is thus physically interesting. It is moreover technically convenient since the formulae simplify 
considerably and a fully analytical treatment can be carried out. 

Using the series expansion for P, cf. Eqs. (13) and (14), we can write the leading contributions to P as 



c+Pi + — P_i + . . . , 



(32) 



where Pi, P_i do not depend explicitly on c+ and are given in Eqs. (15) and (16). In the following, we ignore all 
terms suppressed by higher powers of c+ . 

Inserting Eq. (32) into the expression for the dilaton field, we find at leading order in an inverse c+ expansion 



4^ ^ 2e^*° sinh^(2p) 
4 P^P{ 



Pi 



Pi 



+ ... 



where the dots indicate terms suppressed by higher powers of c+. Defining 



Q 



P-i P'-i 



Pi 



Pi 



and using the identity (which can be easily derived using Eq. (15)) 

2 sinh^(2p) _ 3 
^Pi " 8 ' 

we can write the D3-brane potential of Eq. (31) as follows: 

3A/-2r3(8puv- l)e-«''"^/3 



V 



16 cl 



iV2 



8 c 



(33) 



(34) 



(35) 



(36) 



This potential clearly vanishes in the limit of 



Indeed, in this limit the D3 branc becomes a BPS state 



(thus experiences no force). The function M in Eq. (34), characterizes the field dependence of the potential. Its 
features are controlled by the functions P_i and Q, quantities that control how much the geometry deviates from KS. 
Hence, it is clear that the features of the potential directly depend on the characteristics of the geometry. Using the 
following identity 



P!_i _ 2P_i g2 _^ 4 sinh2 2/5 

'p[~~^^p!^ 



P{P2 



dp 



PjQQ' 
sinh^ 2p 



which can be derived using Eq. (16), we easily obtain an integral expression for A4{p): 



N^M{p) = 



4 sinh 2/9 



\. PjQQ' 



dp — 



PiQQ' 



Pi Pi Jp sinh" 2p 4 
This will be of great use below. In particular, the derivative of M reads 

2 dQ^ 



sinh 2p 



dp 



P2 dp 



NlM\p) = - 



p2 dp 



(37) 



(38) 



(39) 



implying that the force acting on the brane is always attractive since the derivative of as obtained from Eq. (2) is 
always positive. 

At leading order in a c^^-expansion, the relation between the canonically normalized radial coordinate r as given 
by Eqs. (27), (31) and the usual coordinate p does not depend on c+, and reads 



dr ■ 



3 (8puv - 1) 



= ^/2N^ 



2 

2 (8puv 



1) 



VPie- 

sinh 2p 
Pi 



3 

"f"''' dp. 



(40) 



e 
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The formulae become even simpler when focussing on the region p p^. This is the region on which we would 
like to concentrate our attention: precisely in this region the attractive force acting on the brane has a maximum at 
p ^ Pq, and then starts to slowly decrease in size towards larger values of p. 

In order to study this feature, we focus on a limit in which (or more precisely e''* ) is much larger than one. In 
this limit, we can write the leading terms in Pi and its derivative as 



Pi ~ e 



Ap,/3 



1 + e-^P' (sinh(4p) - 4p) 
o 



16 



Thus, we have 



and 



Pic.^e-'f"/'smh\2p), 



N^M'{p) ~ -2e-^P' 



dp ' 



(41) 

(42) 

(43) 



~ 2 



3(8puv-l) 



g-2pw/3 g-4p,/3 sinh(2p) d/9 . 



(44) 



At this point, we have all the tools necessary to analyze the maximum for the attractive force acting on the brane. 
This corresponds to a point in which the second derivative of the brane potential V along the canonically normalized 
radial coordinate r vanishes: it is an inflection point of the brane potential. 

A straightforward calculation shows that such a point exists, and is (approximately) determined by the root of the 
following equation in the coordinate p: 



-5 + 40p^ + 4(1 + 6p^) cosh4p + coshSp - p (30sinh4/9 + sinhSp) = . 



(45) 



Remarkably, this equation does not involve any of the parameters defining the geometry and is valid in the limit in 
which both 0+ and e^* are much larger than one. Numerically, one finds that the root of the previous equation is at 
Pq ~ 0.98 — we have checked using numerical solutions that large but finite values of c+ and p* only lead to tiny shifts 
in the position of po- For p> po, the attractive force acting on the brane starts to decrease in size. In the following, 
we will focus on the region po < P ^ P* in our search for an inflationary trajectory. The effective potential depends 
on the tunable parameters c+, p*, puv that allow us to change both its size and shape. Geometrically, they control 
the profiles for the dilaton, warp factor and the p-form fields. These parameters will provide considerable flexibility 
for tuning the properties of our inflationary model. A qualitative understanding of the shape of the force acting on 
the probe D3-brane using aspects of the dual QFT is postponed to Appendix B 1. 



B. Inflatioriciry quantities and their scaling behavior 

In the limit of large c+ and p* , we can obtain analytic expressions in the parameters controlling the inflationary 

dynamics. These analytic expressions are simple enough to enjoy a clean physical interpretation. 

• Slow-roll parameters: Using the expression for the D3 brane potential in Eq. (36), it is straightforward to 
calculate the slow-roll parameters. At leading order, the large-C-|- expansion of the e-parameter in Eq. (28) reads 



64 V2 n A^c4V8/0uv - 1 VAsinh(2p) dp 
This quantity can be further expanded at leading order in p, as 



(46) 



12 



2 



64 T3 N,cX^8pijv - 1 Vsinh(2p) dp 
V3 Ml iV3 e4puv/3 g-8p,/3 /(2pcoth(2p) - l)(sinh(4p) - 4p) 



(47) 



16\/2 73 4V8puv - 1 V sinh^(2p) 

This expression exhibits a scaling behavior: all the explicit dependence on p is contained in the function inside the 
parenthesis, which does not depend on the remaining parameters that instead all appear in the overall prefactor. The 
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maximum of the function of p inside the parenthesis occurs at po — 0.98. This is precisely the position at which the 
force acting on the brane has its maximum, corresponding to the inflection point for the brane potential — see the 
discussion around Eq. (45). Notice that the value of e is very suppressed by powers of c+ and also has an exponential 
suppression in p*: these small factors make the parameter e very small even near the maximum of the force, around 

For the parameter we have, in the simultaneous limit of large c_|_ and p*. 



Ml 



V3_ 
8n/2 T3 
V3 



p4puv/3 



4%/2 



Ml 
T3 
1 



TV g4puv/3 



1 



92 Q2 



1 



a(sinh(2p))-i dQ' 



T Vsinh2(2p) dp^ sinh(2p) 



dp 



dp 



(48) 



sinh^(2p) 



40p^ + 4(1 + 6p^) cosh4p + cosh 8/9 - p{30 sinhAp + sinh 



which also exhibits a scaling behavior: The expression in brackets depends uniquely on p and is multiplied by a 
coefficient depending on the parameters of the model. Eq. (48) implies that also 77 is suppressed by c-f, but with a 
smaller power than the suppression of e. The expression in brackets vanishes for po ~ 0.98. However, it generically 
becomes large in the region p < po, so in order to have slow-roll inflation we have to focus on the region po < p <ti p^,, 
in which the rj parameter is negative. The latter is a nice extra since observations to favour a red tilt of the power 
spectrum [10]. 

We show a comparison of the above approximate analytical expressions with numerical solutions to the master 
equation in Fig. 6. In this way, one can verify the scaling behaviour and determine the quality of the above approxi- 
mations. As one can see from Fig. 6, the approximation for e is valid for p ^ p*, while the approximation for is a 
good one for all values of p already for p* > 2. 



Mp- 





FIG. 6: Plots of e and rj rescaled by the indicated factors of -j^j^ and -^^j^, respectively, to remove the c+-dependence. 

The curves with different colors correspond to solution obtained by numerically solving the master equation for three different 
values of p, ~ 1 (blue), 2 (red) and 3 (green). The solid, dashed and dotted black lines in the left figure display the approximation 
to e given in Eq. (47). In the right figure, we show the numerical solution for rj for p* ~ 1 (blue) and 2 (red); the curve with 
p, ~ 3 would give the same result as the one with p* ~ 2. The solid black line shows the approximation for rj given in Eq. (48). 



For p < Po, our formulae for the inflaton potential and canonical normalization of the inflaton field are no longer a 
good approximation: the inflaton starts to move faster since the slow-roll approximation breaks down, and the DBI 
nature of the D3-brane action manifests. 

Notice that while the profile for e presents a single hump, the profile for rj in Fig. 6 has a richer behavior; it 
decreases vanishing at po ~ 0.98, then has a minimum at slightly larger value of the radial coordinate (p ~ 1.26), 
and then starts to increase towards a vanishing value for large p. Due to the aforementioned scaling behavior, one 
cannot change the radial profiles of the slow-roll parameters in the present treatment, but only their amplitude. It is 
nevertheless possible to add other sources to the geometry in the form of 'flavour' D5-branes [40], which can further 
modify these profiles (we will briefly comment on this possibility in Section V). 

• Number of e-foldings: The number of e-foldings, at leading order in a slow-roll expansion, can be expressed as 

J Ml J drV \dp) Mp J V2^\dpJ 

4V2 T3 cj V8puv - 1 e-^P"^/3 fP-^f sinh^(2p)rfp ^ ' 

~1MmI (2pcoth(2p)-l)(sinh(4p)-4p) 
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where pinf (lying in the interval po < Pinf <C p*) can be adjusted to provide 60 e-foldings of inflation. Within this 
interval, we are ensured that our formula is correct since we are within the slow-roll regime. Note also that again 
there is a scaling behaviour — the model parameters enter only into the overall prefactor. 

Once properly normalized, pinf is associated to the position of the D3 brane 60 e-folds before the end of slow-roll 
inflation. Comparing the expression for the 77-parameter, Eq. (48), with the previous formula Eq. (49), we find a 
remarkable relation: 

H^^JiP^nf), (50) 

where F{pinf) is a mild function of Pm/, whose values are between zero and (about) one half: see Fig. 7. It is obtained 
by multiplying the absolute values of the two p-dependent functions in Eqs. (48) and (49) since the prefactors depending 
on the parameters in 77 and A^e are precisely the inverse of each other. Hence, for values of pm/ larger than say 3, the 
parameter 77 is ^ — l/(27Ve). 
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FIG. 7: The function F{p) appearing in Eq. (50). 




• Fixing the normalization of the power-spectrum: Besides realizing at least 60 e-folds of inflation, we also 
have to ensure that the amplitude of the scalar power spectrum of primordial fluctuations produced during inflation 
matches the observed value [10]. This amounts to require that the following condition is satisfied 



V2 



V 



Ml 



2.4 X 10" 



V 



2 Mi 



5.6 X 10" 



(51) 



P=P„^f 



By tuning the available parameters of the model (for example the quantity puv) it is not hard to satisfy this condition. 
As we will comment below, the tilt Ug of the scalar power spectrum can also be easily accommodated to fit observations. 

To conclude this section, let us emphasize once again that we have enough free parameters available for tuning the 
size of the slow-roll parameters independently. The geometrical features of the throats we consider are rich enough to 
allow for a flexible implementation of an inflationary scenario. The parameter c+ controls the profile of the dilaton, 
p^ is the position at which the F5 and H3 fluxes get turned on (see also [35]), while the parameter puv controls the 
position at which our throat is glued to a compact space. In the QFT dual, the dynamics in the slow-roll inflationary 
region is mostly controlled by the VEVs of two operators, which are absent in the KS case: a dimension-2 VEV that 
brings the background on the baryonic branch [36] and the VEV for a dimension-6 operator specifically associated 
with the walking dynamics. Although we have freedom to choose these parameters, our scenario leads to robust 
predictions: the slow-roll parameter e is much smaller than the ?7-parameter; the ?7-parameter is negative, and its 
absolute value is proportional to the inverse of the number of e-folds. 



C. An explicit realization of an inflationary trajectory 

Let us apply these formulae to determine a concrete example of an inflationary trajectory in the range po < p ^ p*. 
We focus on p ^ p* since in this region our analytic expressions are accurate. When approaching p ~ pn, our 
approximations break down; in this region the geometry approaches very well the KS throat where the force on the 
brane vanishes. In the region p < poj on the contrary, the attractive force on the brane is typically too strong to 
support slow-roll inflation. The ry-parameter becomes large and the inflationary evolution cannot gain many e-foldings 
in this region. We expect that once the brane enters the region p < po it rapidly falls towards the tip of the throat (at 
p = 0), where it starts backreacting on the geometry making inflation end (see Appendix B 2 for a more quantitative 
treatment of this backreaction) . Hence, we consider p ~ po as the lower value of the radial coordinate at which 
standard slow-roll inflation ends. 
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Given these considerations, we would like to determine an inflationary trajectory generating sixty e-folds of slow-roll 
inflation nearby the inflection point for the potential, namely for p ^ po- As we learned in the previous section, the 
slow-roll parameters are suppressed by powers of Yc+ and e~''* . In a regime of small values for these quantities, 
as the one we are interested in, it is not difficult to adjust the slow-roll parameters 77, e to sufficiently small values 
and moreover tune them independently, maintaining the robust relation e <^ \ri\. The normalization of the power 
spectrum, then can be easily further tuned by choosing puv- (Note that the scalings of all quantities depend only on 
certain ratios of c-|_ and Nc such that any change in can be compensated by adjusting c+ accordingly.) 

Before providing the actual numbers, let us give an argument to fix the value of T3, the tension of the moving D3 
brane, which depends on the details of bulk moduli stabilization. Upon introducing a finite UV cutoff puv and gluing 
the throat into a compact space, the graviton zero mode will have most of its support in the essentially unwarped 
compact space since it is exponentially suppressed in strongly-warped regions (see e. g. the discussion around Eq. 
(C.IO) of [3]). That is, we can approximate the (four-dimensional Planck-mass) Mp as 



Ml 



2Ve 



(52) 



where Vg is the volume of the compact space. 

Recall that when defining the backgrounds around Eq. (1) we have set a'gs = 1. This was used in deriving the 
above expressions for N^, e and 77, so now we have to reinstate appropriate powers of a'g^ to make iVe, e and 77 
dimensionless. One can easily see that the relevant combination that should appear is 



T3 a'gs 
Ml 



(53) 



where the D3-brane tension Ta is given by 



(54) 



(27r)3a'2g, ' 

and Mp is determined by Eq. (52). For the numerical estimates in the following, we assume the bulk moduli controlling 
gs and Vg have been stabilized such that (/s ~ 0.1 and Vg ^ ^5 \/a'^ , which yields 

T^a'g, {2^ya'^gl a'g, (2^)^ 



Ml 
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(55) 



Similarly, the ratio T3/MI is given by 



Ml 



(27r)^a'6 
AVI 



6 X 10" 



(56) 



Equipped with these two ratios (the latter one enters into the amplitude of the scalar power spectrum) , whose precise 
values depend on the details of bulk moduli stabilization, we can now fix our parameters in order to match observations. 






FIG. 8: Plots of e, rj and Ne for the example choice of parameters specified in Eq. (57). The gray dashed lines mark the 
positions of the inflection points po — 0.98 and the point pi^f corresponding to about 60 e-folds before the end of inflation. 
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Wc now choose the remaining parameters in order to reahze inflation — with mild tunings one easily finds a 
satisfactory inflationary trajectory. For example, the following choice works: 

p* 2^2.90, /Ouv=9, c+~ 6.1x103, iVc = 1 , p™/ = 1.8 (57) 

which yields Ne ~ 61.5 and matches the observed power spectrum normalization. At the value pinj = 1.8, that is 
sixty e-folds before the end of slow-roll inflation, we find e ~ 7 x 10"^^ and rj ~ —6 x 10^"^. Thus, the value for the 
scalar spectral index ng = 1 + 2r7 — 6e ~ 0.988 which is consistent with the value observed by WMAP within 2ct [10]. 
The tensor-to-scalar ratio r = 16e ^ 10""'^^ is completely negligible in accord with the non-observation of primordial 
gravitational waves; this is a generic feature of warped D3-brane inflation, cf. e. g. [16-19]. See Fig. 8 for plots of e, 
T] and A^e as a function of pj„/ for the above choice of model parameters. 

As we explained above, the ry-parameter increases and becomes large for p < po! this implies that one does not gain 
many e-foldings while moving in this region. On the contrary, the e-parameter stays very small, and this would seem 
to imply that inflation cannot end in this set-up. However, as we discuss in Appendix B 2, when the brane reaches the 
deep IR, its backreaction on the geometry becomes important and our description of the throat is no longer reliable. 
In this regime deformations of the space-time should be so drastic to terminate inflation. 

V. DISCUSSION 

We have realized a model of slow-roll D-brane inflation in a warped geometry dual to a fleld theory exhibiting 
walking behavior: hence its name slow-walking inflation. In our set-up, the rich fleld content of the geometry we 
consider is associated with various parameters that control an analytically calculable inflationary potential. It has 
the correct properties to drive sixty e-folds of inflation and to match the observed values for the amplitude and scale 
dependence of the power spectrum of scalar perturbations. Interestingly, the D-brane potential generically has an 
inflection point, induced by the particular shape of the dilaton profile. The relevant features of the geometry have 
a counterpart in a dual field theory, allowing to understand the characteristics of the inflationary brane potential in 
terms of the combined action of operators with different dimensions. Our model can be seen as a generalization of 
the inflationary set-up in the baryonic branch of the KS throat [9]: with that scenario, it shares the good feature of 
not needing an anti-D3-brane to drive inflation, since the warping and the fluxes turned on act with sufficient force 
on the brane to drive inflation. Moreover, inflation naturally occurs in a region that is easy to identify from the 
geometrical point of view, far away from the place in which the throat is glued to a compact space (as needed to 
couple gravity to the system and to add ingredients for stabilizing light moduli on the throat). Our methods can be 
straightforwardly applied to warped geometries that generalize the ones we considered here, containing e. g. additional 
source D5-branes with non-trivial proflles in the throat [40-42]. This allows to further enrich the potential for the 
moving probe D3-brane, providing additional freedom to tune and engineer the inflationary process. It would also be 
very interesting to study dynamics of the brane along angular directions [43] , and analyze whether our set-up can be 
used to realize DBI inflation in the region in which the //-parameter becomes large [44]. We leave these questions for 
future work. 

Let us make some final caiitionary remarks about the 77-problem. The moving D3 brane couples with the background 
fields, hence it backreacts on the geometry. In Appendix B 2 we show that the backreaction on the throat fields is 
mild, until the brane reaches the far IR region. However, once the compact CY is glued to the throat, new couplings 
between bulk objects (introduced to stabilize light moduli) and the moving D3 brane itself, can modify the brane 
potential. All this can lead to order one corrections to the infiationary 77-parameter. The effects of gluing the compact 
CY to the throat can be seen also from the point of view of the dual field theory: new relevant operators become 
important in this case. For a complete model, it will be necessary to study the properties of these contributions. This 
is not easy to do and we leave it as an open problem for the future. Techniques as the ones developed in [13-15] 
should be used for this purpose. Nonetheless, as we discussed above, our model has some new promising attractive 
features worth to study further. It exploits the rich geometrical content of the throat we considered to build a model 
of inflation, in which the brane potential is fully calculable. The resulting features of the inflationary potential can 
be tuned with much freedom thanks to the many parameters controlling the; geometry. 

Optimistically, this might imply that even after including the aforementioned corrections due to physics in the 
compact space, there is still room left for changing the parameters characterizing our original model in such a way 
that sufficient slow-roll inflation can be achieved without qualitatively changing much the framework we developed. 
We leave these interesting questions for future work. 
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Appendix A: Details of the backgrounds 

In this Appendix, we give some details about the backgrounds sketched in Eqs. (1) and (6). As explained in the 
main text, we consider the geometry produced by stacking on top of each other Nc D5-branes that wrap an S'^ inside 
a CY-cone, e. g. the conifold. Let us define a set of SU (2) left- invariant one-forms as 



(Di = cos ipd6 + sin ^ sin Odcp , 



UJ2 = — sin tpdO + cos tp sin 6di^ , 



u>3 = dip + cos 6di^ . 



(A.1) 



The conifold with topology R x S'^ x has a SU{2) x SU(2) x C/(l) isometry. 

We assume that the functions appearing in the background depend only the radial coordinate p (the range of the 
angles isO<9,9<Tr,0<(f>,(f)<2Tr,0<tp< 47r). We write the background (in Einstein frame) as. 



ds' 
dsl 



e*''')/^ [dxl. 



d4] , 



„2fc 



»2«(p) 



{d0^ + sin^ 0d(t>'^) 



+ — — [(wi + a{p)def + {uj2 - a{p) smOdcpy] + — j— (^s + cos 9d^y , 



4 



(A.2) 



(wi -I- b{p)d9) A {Cj2 - b{p) sin9d(j)) A (0)3 -|- cos 



+ dpb dp A {-d9 Aui + sin 9d(p A W2) (1 - b{p)^) sin 9d9 A A W3 



The full solution is determined by solving the equations of motion for the functions {a, b, g, q, k}. As described in 

the main text, the BPS equations derived using this ansatz can be rearranged in a more convenient form, by rewriting 
the functions appearing in Eq. (A.2) in terms of a new basis of functions where the BPS equations decouple, as 
explained in [32], in which the functions {a, b, g, q, k} read 



4e2« = „ ^^^'^^ ^ , e^s = P coth(2p) - Q , 



„2k 



Pcoth(2p) -Q 
P' 

Y' ' 



sinh(2p) (Pcoth(2/5) - Q) ' 



2p 



(A.3) 



sinh(2/9) 



Using these new variables, one can work with the BPS equations to obtain a single second order equation for P{p), 
while the quantities Q and $ are given by 



Q{p)^NA2pcoth{2p)-l), 
4^ _ 8e*'^o sinh^(2p) 

with $0 a constant and P determined by the second order equation mentioned above: 

P" + P'r^ + ^-4coth(2p))=0. 



(A.4) 



(A.5) 



In the case of the background after the rotation, the details go as follows. We use the vielbein in Einstein frame given 
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by Eq. (4). The full generated configuration is obtained as 

10 



i=l 

_ 3 

e 4 



hi 

3*/2 

B2 = k2^— [e"^ + cosM (e^^ + e'^) + sin/x (e^^ + e^^)] , 



5*^ 

e 4 



'dp 



where 



P coth(2p) - Q 
the functions fi,i = 1 , . . , 4 are given by 

/i = -27V,e-'=-2s , f^ = ^{a^- 2a6+ l)e-'=-23 ^ 



(A.6) 



P-Qcoth(2p) 



(A.8) 



and we denoted 

eijk...i = A e^' A A • • • A e' . (A.9) 

A necessary condition to apply this solution-generating technique is that the quantity e* is bounded from above, being 
an increasing function with e*(°°' its maximum value. This condition can be linked with the absence of D7-brane 
sources in the configuration of Eq. (A.6) — see [34] for details. 



Appendix B: Some comments on the dual QFT & brane backreaction 

Here we collect some arguments that should be thought as 'plausibility arguments' for the results obtained in this 
paper. The first of these arguments, which is of qualitative nature, attempts a QFT explanation for the form of the 
potential and the force felt by the travelling D3-brane. The second argument aims at justifying that our set-up, with 
an inflection point around which inflation starts, does not get substantially modified when including in the effects of 
back reaction of the travelling D3-brane on the geometry. 



1. A fleld-theoretical view on the form of the potential and the force 

The goal of this Appendix is to get some handle on the force acting on the moving D3 brane using arguments in 
the dual QFT. The argument we will present is of qualitative nature. We will basically only attempt to make the 
point that the 'shape' of the force on the D3 as computed in the gravitational set-up has some connection with the 
same quantity calculated in a QFT that approximately describes the process. 

In order to make this discussion more self-contained, we will summarise some aspects of the minimally supersym- 
metric QFT associated with the rotation of the walking solution as introducted in Section II B. 

Let us start by considering the quiver field theory 

SU{{k + l)Nc + nf) X SU{kNc + n/) , (B.l) 

where in this case n/ = 1 represents the single D3- brane that is moving along the throat, representing our candidate 
inflaton field. 



19 



For the case Uf = 1, the field theory was weh-studied by Seiberg in the 1990's (see for instance [45]). In our case, 
we consider a supersymmetric QFT with quiver 

SU{{k + l)iVc) X SU{kNc) , (B.2) 

and add one probe D3-branc. This breaks SUSY generating a non-constant potential and a force on the probe D3. 

This type of SUSY breaking may be thought of as being 'soft'. We will therefore use Seiberg's supersymmetric 
formalism, knowing that the modification due to soft breaking will not change the qualitative results — in other 
words, in our model, we assume that the scale of SUSY breaking is smaller than the strong coupling scale of the QFT 
in Eq. (B.2). This was done successfully in this context in the papers [46]. 

Let us go to the point now and summarise some things about the QFT's described above, adapting arguments of 
[9]. The correct indices in the following expressions are explicitly written in [9]. 

• The field theory of Eq. (B.l) has a superpotential given by a tree- level part and a non-perturbative part of the 
form 

Wr~.M^ + j {MAB - det[M]) , (B.3) 

where the first term corresponds to the tree-level superpotential 

Wiree ABAB - MM . (B.4) 

Here, A, B are the usual bifundamentals joining the two gauge groups and M is the meson superfield dual to 
the position of the D3-brane probe [9] . A, B are the baryon and anti-baryon fields, that in the case of the quiver 
in Eq. (B.2) are given by 

A = e''-^''Ai,...Ai^, B = e''-'''Bi,...Bi^. (B.5) 

For the quiver in Eq. (B.2), N = 2Nc at the end of the cascade where the quiver becomes SU{2Nc) x SU{Nc) 
with the SU{2Nc) factor being strongly-coupled. In the case of the quiver in Eq. (B.l), the last step of the 
cascade gives the gauge group SU(2Nc + rif) x SU{Nc + Uf). In the strongly-coupled sector, the number of 
colors is 2Nc + n/, while the number of fiavors is 2Nc + 2nf. Since the index N in Eq. (B.5) runs over fiavor 
indices, then for n/ = 1 (one probe D3-brane), the baryon operators have one free fiavor index. 

• For simplicity, we assume that the Kahler potential is trivial. This implies that up to numerical factors, the 
potential reads 




(B.6) 



• We will also consider that the difference between the VEVs of the baryon and anti-baryon fields is related to 
the VEV of the operator U, defined in [9], 

Ur^{A-B). (B.7) 
This is the case for example in the KS-limit of the quiver, where U = 0. 

• A crucial point is that the operator U is related to an invariant of the rotation procedure described in Section II B. 
This invariant is the function 

F = a^ + 4e2''-2s _ i , (B.8) 

where the functions a{p), h{p),g{p) arc defined in Eq. (A. 3). The relation between U and this particular com- 
bination of gravity fields above is strengthened by the fact that when analysing the asymptotics of the field F 
above, it is found that it describes a dimension-2 operator getting a VEV [35]. The function F, an invariant 
under the rotation procedure, is drawn in Fig. 3 of [35]. 

• We also assume that the meson matrix has one non-zero component while the remaining components vanish, 
such that det M = 0. 



20 



Under these assumptions, we now calculate the potential from Eq. (B.6) and obtain 

V ^ {M - jABf + M^A"^ + B^) . (B.9) 

From this we may calculate the force on the meson field /m ~ —OmV. We consider the situation in which (M) 
changes slowly on a background of fixed {A) , {B) , we then have 

fw- A^ ^B"^ -AB^ {A-Bf ^U'^ (B.IO) 

Let us summarise: we are proposing a very qualitative reasoning to get a handle on one of the scales of the problem 
and its effect on the force on the probe D3. Prom the QFT perspective, we are dealing with a case of very soft SUSY 
breaking. For a field theory like the one in Eq. (B.l), for nj = 1, we have a potential like in Eq. (B.3). But we are 
studying this by starting from a theory like in Eq. (B.2) and using the fact that this theory has VEV's for the baryon 
and anti-baryon superfields. Then we break SUSY 'softly' by the addition of a moving D3 brane. Hence, wc propose 
that wc have 'background' baryonic fields taking a large VEV, while the mcsonic supcrficld probes this QFT, varying 
slowly. The result of this very qualitative reasoning tells us that the force on the D3 probe has a similar shape as 
that of the 'geometry' invariant given by the combination in Eq. (B.8). In terms of the radial coordinate, this is a 
function which has a maximum close to an intermediate point — that turns out to be — and then decays in the 
IR and in the far UV. Notice that the QFT dual to our geometry has, as discussed in [35], two independent scales. 
One of the scales controlled by indicates, when flowing from the UV to the IR, the beginning of the walking region. 
Flowing further towards the IR, the second scale is given by p/ ~ 1. This scale indicates the end of the walking region, 
the onset of confinement and 'L2N^ symmetry breaking. These two scales, not surprisingly control our inflationary 
dynamics. 

Note that we obtain a qualitatively very similar result for the e slow-rolling parameter in Eq. (28). This is due to 
the fact that roughly e ^ {¥' jV)'^ , for an almost constant V . In this sense, there should be some connection between 
the QFT potential and the potential felt by the D3 probe. Note that the point in the radial direction where the 
force is strongest corresponds to the transition between the QFT being dominated by its UV behaviour to it being 
dominated by its IR behaviour. This transition point dominates the behaviour of our model of inflation. 

Let us move now to a more quantitative and precise way of understanding the effects of the backreaction of the D3 
brane on the throat. 



2. Backreaction of the travelling D3-brane 

The aim of this subsection is to explore how good the approximation of a travelling ■prohe D3-brane, actually is. 
The outcome of the study is that, for slowly moving D3 branes, the backreaction effects are important only close 
to the far IR region p ~ 0. In other words, we will give a quantitative way of calculating the backreaction of the 
D3-brane on the geometry, in a particular situation — this will be in concordance with the approximations we used 
in our model of inflation. 

To make this calculation feasible, wc assume for simplicity to be in a KS-background (since we are taking very 
large values of c+ this is not a very drastic assumption) and add D3 sources. We will also assume that the D3-brane 
moves quite slowly. That is, its kinetic energy is negligible when calculating the deformation of the space time that 
is mostly sourced by its mass. 

The main simplification is to consider the backreaction assuming the D3 to be smeared on the directions 0, </?, 0, (p, ip. 
If we focus on this particular and simplifled problem, we can write an exact solution. The solution is actually written 
in Eqs. (3.8) and (3.9) of [40]. There, it was found that the warp factor h becomes 

i, . , 4fci Sjx) ^ 25/3fe^^2 .CO (sinh(4a:) - 4a:)V^(2xcoth(2a:) - 1) 

where S{x) is the probe D3 distribution. We should take S{p) = 6{p — po) or some other function localized close to 
the position of the probe D3, po, say like a gaussian or similar. The point po (which should not be confused with the 
position of the inflection point introduced in Section IV A) slowly changes as the brane moves towards the end of the 

geometry. 

Suppose that we want to estimate the backreaction for large values of the radial coordinate. In this case the warp 
factor will change as 

h-fiKS oc r • u^^ — 7 — A — TTT^ JTT^ (B-12) 

[smh(4po) - 4po] ' e'^P°'^ 
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where we approximated the profile of branes by a delta function localized at pQ. Thus, the correction is very small. 
However, for the case in which the brane is close to p ~ 0, we have that the new contribution to the warp factor goes 
like 

h- fiKS y^-tTa — r — A — iTrt ~ ~~ (B.13) 
[smh(4/>o) - 4po]^/^ Po 

for small po- Hence the correction is large. 

Now, to the point: for values where inflation occurs, we are able to use the approximation above and see explicitly 
that the D3 probe should actually not change things too much. After the end of Inflation, we come close to p = 
and the set-up and the solution should change due to the backreaction. 

Things should not be very different for large but finite values of c+, in which case an exact solution cannot be 
written, but the system is close to KS. 

There is a field theoretical way of stating the argument above: if the meson field is 'slowly varying', that is the D3 
is slowly rolling, we can integrate out the meson in the super potential of Eq. (B.3) and we would then get the super 
potential that one typically has in a quiver of the form in Eq. (B.2), giving the usual baryonic branch as one of its 
vacua. The gravity solution is not much different to the (baryonic branch of the) Klebanov-Strassler background, this 
is expressed by Eq. (B.12). As the meson VEV becomes comparable to the other scales in the problem, the meson 
field becomes more influential and we should instead describe the SUSY dynamics of the quiver in Eq. (B.l). This 
coincides with the D3-brane going to the far IR of the geometry and changing it considerably — this is expressed by 
Eq. (B.13). In all the equations above, we were assuming a SUSY situation (or very 'softly' broken SUSY as in [46]). 
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